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Experimental realization of Josephson junctions for an Atom SQUID 
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We report the creation of ideal Josephson junctions in a toroidal dilute gas Bose-Einstein conden- 
sate (BEC). The demonstrated configuration of a pair of junctions on a multiply-connected BEC is 
the cold atom analog of the well-known dc Superconducting Quantum Inteference Device (SQUID). 
We measure the critical current of the junctions, observe Josephson effects, and find dynamic behav- 
ior that is in good agreement with the simple Josephson equations for an ideal tunnel junction with 
a sinusoidal current-phase relation. The junctions and toroidal trap are created with the Painted 
Potential, a time-averaged optical dipole potential technique which will allow scaling to more com- 
plex BEC circuit geometries than the single Atom-SQUID case reported here. Since rotation plays 
the same role in the Atom SQUID as magnetic field does in the dc SQUID magnetometer, the 
device has potential as a compact rotation sensor. It may also be useful for creating macroscopic 
Schrodinger Cat states. 



The dc Superconducting Quantum Interference de- 
vice (SQUID) 1 was created to study and utilize 
the quantum interference of currents flowing in parallel 
through Josephson junctions connected in a supercon- 
ducting loop. Today, dc SQUIDs, and the more recently- 
developed rf SQUIDs, are the basis of some of the most 
sensitive magnetometers [T]. An "Atom SQUID" is an 
analogous quantum interference device that uses a su- 
perfluid dilute gas Bose-Einstein condensate (BEC) [21 E] 
flowing through potential barriers in a trap. The su- 
perfluid counterpart of the SQUID response to magnetic 
fields is an analogous response to rotation, which has 
already been demonstrated in the laboratory with su- 
perfluid helium analogs of the SQUID having sensed the 
rotation of the Earth [IHS]. The Atom SQUID is partic- 
ularly attractive for rotation sensing and for addressing 
basic questions in quantum physics because of the advan- 
tages of easy detection of BEC atom number and phase, 
along with the existence of accurate well developed mi- 
croscopic theories of condensate physics. 

An Atom SQUID has two essential ingredients: a 
multiply-connected geometry trap with quantization of 
superfluid circulation, and potential barriers (junctions) 
exhibiting Josephson effects. Like SQUIDs, Atom 
SQUIDs can have two different types: a "dc Atom 
SQUID" with two junctions, and an "rf Atom SQUID" 
with a single junction. Previously, we created a BEC in 
a toroidal trap and demonstrated quantization of circu- 
lation in that system [7] . Here we report the first experi- 
mental creation of Josephson junctions in a toroidal trap 
and the observation of Josephson effects with them. Both 
the double Josephson junctions and the toroidal trap for 
a dc Atom SQUID were created with the "Painted Poten- 
tial" technique [5] . This technique is general and may be 
used to create arrays of Josephson junctions with com- 
plex spatial and temporal variations. To demonstrate 
Josephson effects, a bias current was imposed by moving 
the two junctions relative to each other. Compression of 
atoms was observed when the bias current exceeded the 
critical current of the junctions, allowing the critical cur- 



rents to be measured. This advance opens the door to the 
study of Atom SQUID physics, including the investiga- 
tion of quantum interference of currents due to rotation 
and its application to rotation sensing. 

Creating an Atom SQUID is a challenge because both 
the toroidal trap and the Josephson junctions are com- 
plex potentials that go beyond conventional magnetic or 
optical dipole trapping geometries. BECs have been cre- 
ated in toroidal traps using a magnetic trap with a repul- 
sive optical trap at the centre [5], a Laguerre-Gaussian 
mode of a laser beam ^U\, |TT] , and a Painted Potential 
[21 H]. Separately, Josephson junctions for BECs have 
been realized with interference of laser beams |12j , a mag- 
netic trap with a repulsive barrier from a laser beam |13| , 
and a radio frequency dressed magnetic trap [13]. Re- 
cently [ini [H] , a separate laser beam was used to create 
a barrier for atoms in a toroidal trap, realizing an atom 
circuit for studying the interesting physics of phase slips 
[TB] in a BEC. That toroidal atom circuit with a single 
rotating barrier is an atomic analog of a "single junction 
dc SQUID", a novel configuration which has not been 
studied in detail in the context of superconductors. The 
thick (compared to the healing length) and low (50% of 
the chemical potential) barrier used in that experiment 
did not form an ideal Josephson junction, and the ex- 
periment did not directly demonstrate Josephson effects. 
In this paper, we report the first experiments showing 
Josephson effects for a BEC in a multiply-connected ge- 
ometry. The junctions in the toroidal trap, realized with 
the Painted Potential technique, are thin and high po- 
tential barriers with significant quantum tunneling. As 
a result they behave as ideal Josephson junctions with a 
sinusoidal current-phase relation. 

The Painted Potential technique [3 |H] realizes com- 
plex dynamic potentials for BECs by rapidly moving 
far-detuned laser beams to create time-averaged optical 
dipole potentials. Improvements to the setup described 
in Ref. 8J, in particular a new long distance microscope 
objective (N.A. = 0.4), make the spatial resolution of 
the Painted Potential (1.5 /im) comparable to the heal- 
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FIG. 1: (a) Schematic of the experimental set up. Arbitrary 
and dynamic potentials for the BEG are created by the super- 
position of two red-detuned beams making attractive optical 
dipole potentials. One beam is a horizontal light sheet, the 
other is a rapidly moving focused beam which paints the po- 
tential, (b) Galculated intensity distribution of 8 ^m diameter 
painted Atom SQUID potential, (c) In situ absorption image 
of a BEG in the Atom SQUID potential. 



ing length of a BEC so that painted potential barriers 
can exhibit significant tunneling rates. Figurema) is 
a schematic drawing of the setup. A horizontal scan- 
ning beam creates a flat two dimensional (2D) potential 
for supporting atoms against gravity, while the vertical 
painting beam creates a complex dynamic potential for 
a BEC. This high resolution system enables painting of 
many variations of complex Josephson junctions, includ- 
ing the dc Atom SQUID configuration of double Joseph- 
son junctions on a toroidal trap. FigureflVb) shows the 
calculated beam intensity distribution for this case where 
the intensity is lowered in two symmetric locations to cre- 
ate tunneling barriers for a BEC, and Fig.lljc) shows the 
in situ absorption image of a BEC in this potential. The 
FWHM (full-width at half-maximum) of the barrier is 
2 /im, which is small enough for a BEC to have a signifi- 
cant tunneling rate (of order a few hundred Hz) . 

It is hard to observe Josephson oscillations in this 
system directly because of their small amplitude. We 
therefore used the scheme suggested in Ref. [l^ to study 
Josephson effects. In our implementation of this idea, 
moving the junctions circumferentially towards each 
other leaves the BEC density unchanged as long as atoms 




FIG. 2; (a) Initial configuration of Atom SQUID potential, 
(b) Final potential after the movement of the two Josephson 
junctions, (c) Image of a BEG where bias current is less 
than critical current, (d) Image of a BEG where bias current 
exceeds critical current, showing the BEG compression which 
occurs when tunneling is not possible. 



can tunnel through the junctions to maintain the same 
chemical potential in both sectors of the torus. This phe- 
nomenon of a current of atoms flowing without a chemi- 
cal potential difference is the dc Josephson effect in this 
system. The current will increase with barrier velocity 
until the critical current of the Josephson junctions is 
reached, at which point the system switches to the ac 
Josephson regime. Here there is an oscillating current of 
atoms through the barrier, with frequency proportional 
to the chemical potential difference across the junctions, 
but no net current across the junctions. Therefore at 
barrier velocities greater than the speed limit imposed by 
the critical current, the moving Josephson junctions sim- 
ply push the atoms, resulting in compression of atoms in 
one side and expansion in the other. Figure^ a) and (b) 
show the initial and final potentials of the Atom SQUID 
for the movement of Josephson junctions. Figurel2[c) 
shows the unchanged density distribution of atoms when 
atoms tunnel through barriers to maintain equilibrium 
and Fig.[2Fd) shows the compressed cloud of atoms in 
one side when the junctions move too fast for tunneling. 
This data is the first demonstration of Josephson effects 
in a toroidal BEC. The rest of the paper discusses the 
relevant theory and the experimental results in detail. 

The simplest theoretical model for the experiment is 
based on the Josephson equations developed for a BEC 
[171 [18]. Although our Atom SQUID setup has two junc- 
tions, because the two junctions are symmetric and there 
is no rotation the system can be modeled as two conden- 
sates connected by a single junction. Letting Ni,N2 and 
01 , (t>2 be the number and phase of atoms in the two sides 



of the torus, and defining the relative population differ- 
ence z = (TVi — N2) / {Ni + N2) and the relative phase 
difference (f) — (f)! — it'2, the Josephson equations are 
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Here Ic — 2Ej /UN is the critical normalized current of 
the double junction, where Ej is the Josephson coupling 
energy [T71 [TO] and A^ = A^i + iV2 is the total number of 
atoms. Also, loc — Ec/2hN, where Ec is the capaci- 
tance energy [T71 [H] , and zq is the population difference 
in the equilibrium condition where the chemical potential 
difference 
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is zero {fii and ^2 are the chemical potentials in the 
two sides of the torus). Equations (1) and (2) determine 
the evolution of the number and phase differences from 
given initial conditions. The derivative i is a normal- 
ized atom current, with maximum value Ic- Changing 
the equilibrium population difference zq creates an effec- 
tive bias current zq. Our experiment follows the proposal 
by Giovanazzi et al. |17| to create Zq through the rela- 
tive movement of Josephson junctions against a BEC. 
This movement changes the trapping volumes in each 
side [Fig.l2ja) and (b)] and hence the number difference 
at equilibrium, which creates the desired bias current. As 
discussed above, zq < Ic is the dc Josephson regime in 
which the chemical potential remains the same in both 
sides of the torus, while for zq > Ic (ac Josephson regime) 
there is no dc current and so the atoms are compressed 
in one side and expanded in the other. 

We now discuss details of the experiment. A BEC 
of ®^Rb atoms was prepared in a potential with trap- 
ping frequencies of 300 Hz in the vertical direction of the 
horizontal trap and 570 Hz along the radial direction of 
the toroidal trap. The preparation of the BEC starts 
with the trapping and evaporative cooling of atoms in 
a quadrupole magnetic trap. ^^Rb atoms in the F = 1 
and mp = — 1 state were cooled to 30 fiK and trans- 
ferred to the combined optical trap. The horizontal trap- 
ping beam, with wavelength A = 1064 nm and waist 
luq = 11 /im, was scanned horizontally at 14 kHz over 
a width of 47/Lim, creating a flat 2D potential for the 
atoms. The vertical trapping beam, with A — 830 nm and 
Wo = 1.5 /im, painted the potential at a scan frequency 
of 25 kHz, creating here an 8 /im diameter ring with two 
intensity minima forming Josephson junctions. Evapora- 
tion in the combined optical trap was accomplished by 
lowering the horizontal trap depth from 40 /iK to 1.1 /iK 
in 2.5 s while the depth of the vertical painting poten- 
tial remained fixed at 75 nK. The result was a pure BEC 



with atom number A'' set between 1000 and 8000, corre- 
sponding to chemical potentials from 29 nK to 58 nK. To 
prevent excitation of flow of atoms in the BEC creation 
process, the initial potential had weaker barriers (22 nK). 
This was sufficient to prevent any rotation of the atoms 
in the toroidal trap throughout the evaporation process. 
At the end of the evaporation stage the barrier height 
was increased to the final value of 44 nK in 100 ms. The 
two junctions were then moved towards each other at a 
rotation frequency / for each junction, which gives rise 
to a bias current zq = 4/. To prevent plasma oscilla- 
tions and to establish the bias current adiabatically, the 
barriers were accelerated at constant rate to the desired 
constant velocity. The total compression angle for each 
junction was tt/S and the compression angle after the ac- 
celeration stage was O.OGtt. Following the movement of 
the Josephson junctions, the BEC was imaged by absorp- 
tion in situ. The clear change seen in the final density 
distribution between the dc Josephson [Fig.[2Fc)] and ac 
Josephson [Fig.[2ld)] regimes as the bias current increases 
can identify the critical current accurately. 

Predicted critical currents were obtained from Joseph- 
son coupling energies computed by applying the im- 
proved two mode model of Ananikian and Bergeman |19j 
to accurate numerical wavefunctions generated through 
direct minimization of the appropriate Gross-Pitaevskii 
energy functional. The results are plotted in Fig. [3] for 
three different potential depths, showing that the criti- 
cal current may be varied substantially over the range of 
atom numbers accessible to the experiment. Figure[3]also 
illustrates the dependence of critical current on potential 
depth, which arises because the critical current depends 
on tunneling (the Josephson coupling energy Ej) which 
in turn depends on the barrier height, which is here a 
fixed fraction of the trap depth. 

Data was collected by fixing the bias current and vary- 
ing the atom number. It is straightforward to determine 
the total number of atoms N and the normalized number 
difference z from the absorption images. Figurep[a)-(c) 
shows three data sets corresponding to barrier rotation 
frequencies of Is"""^, 2s~^, and 4s~^. The sharp tran- 
sition from dc to ac Josephson regimes is seen as the 
sudden change of normalized population difference z as 
the atom number N decreases. The transition point, cor- 
responding to the atom number where the bias current is 
equal to the critical current, shifts to higher N as the bias 
current increases, in accordance with Fig.[3J To make 
a quantitative comparison with theory we find the best 
least-squares fit to each data set of a numerical integra- 
tion of the Josephson equations (1) and (2). Each fit had 
three free parameters: the potential depth, and the initial 
and final values of the normalized population difference 
z. These parameters are determined by the fit because 
slight deviations of the initial and final population differ- 
ences from the values predicted by geometry result from 
imperfections in the toroidal potential, and drifts of a 
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FIG. 3: Solid lines show computed critical current as a func- 
tion of BEG atom number for three different toroidal trap 
depths. Points show the measured critical atom number de- 
termined by fitting a Josephson equation model to data sets 
for the three currents employed in the experiment (see text). 
Critical currents and fitted trap depth are respectively: 1 s~^ 
and 71nK (circle), 2s~-^ and 74.5nK (triangle), 4s~-^ and 
75.6 nK (square). 



few /ini in the relative waist positions of the horizontal 
and vertical trapping beams can introduce smaU varia- 
tions in trap depth between data sets. The fits, shown as 
the blue curves in Fig.|4J are in generally good agreement 
with the data. The potential depths obtained from the 
fits (71 nk, 74.5 nk, and 75.6 nK respectively) are within 
the range expected from our estimate of possible relative 
drift of the waists of the trapping beams. The critical 
atom numbers determined by the fit for each rotation 
frequency are shown on Fig. [3] Their consistency with 
the theory predictions (for the trap depths determined 
by the fits) illustrates the agreement between experiment 
and theory seen in Fig.[4ja)-(c). 

The model based on the Josephson equations (II]) and 
(pi) is simple, neglecting, for example, changes in the 
coupling and capacitive parameters as the barriers move 
and also any consequences of excitations in the conden- 
sate. To quantify the impact of such effects we per- 
formed full three dimensional (3D) dynamic simulations 
for parameters appropriate to our experiment using the 
Gross-Pitaevskii equation (GPE). The results, shown in 
Fig.llFd), are in excellent agreement with the predictions 
of the simple Josephson equations. This result confirms 
the validity of the Josephson equations for describing 
our experiments, and it shows that the Atom SQUID 
can be well understood from microscopic theories with- 
out relying on phenomenological models. This finding 
may be beneficial in the design and interpretation of 
future experiments with Atom SQUIDs. Further, the 
agreement between the GPE results, the data, and the 
Josephson equations confirms that the experiment real- 
izes ideal Josephson tunnel junctions obeying the sinu- 
soidal current-phase relation given by Eq. (fTl). It follows 
that the system reported here is exactly analogous to the 
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FIG. 4: (a) - (c): Measurement of critical currents. Each 
data point is the mean of several measurements of the final 
normalized population difference z, with the error bar being 
the standard error of the mean. The blue curves show the best 
fit numerical integration of the Josephson equations. For the 
three data sets, the junction rotation frequency /, correspond- 
ing bias current i, and duration of junction movement T are: 
(a) / = 0.25Hz, i = ls"\ T = 370ms. (b) / = 0.5Hz, 
i = 2s-\ T = 185ms. (c) / = IHz, i = 4s-\ T = 92.5ms. 
(d) Predicted final value of normalized population difference 
z for / = 1 Hz (bias current i = 4 s~^) in a flat trap computed 
using numerical integration of the Josephson equations (blue 
line) and the 3D Gross-Pitaveskii equation (red disks). 



ideal dc SQUID, and so superfluid analogs of the phe- 
nomena seen in that device should also be observable in 
the device described here. 

In summary, we have demonstrated Josephson effects 
in a toroidal trap for the first time, by showing the tran- 
sition between the dc and ac Josephson regimes and by 
measuring the critical current of the junctions. The ex- 
perimental data for our dc SQUID geometry is in good 
agreement with the predictions of the Josephson equa- 
tions. Looking to the future, the quantization of circula- 
tion in the toroidal trap geometry will lead to quantum 
interference of currents in the dc Atom SQUID as a func- 
tion of rotation rate. The interference can be demon- 
strated by measuring the shift of critical current with 
trap rotation rate (rotating the junctions is straightfor- 
ward with the Painted Potential technique) , enabling the 
use of Atom SQUIDs for rotation sensing. The Atom 
SQUID interaction region can be compact, so it will be 
interesting to compare its rotation-sensing performance 
with the current state of the art, the large area atom 
optics-based gyroscope [20]. Another future direction, 
relevant to both sensing and basic quantum physics, is 
applying the technology demonstrated here to recent pro- 
posals pilES] to create Schrodinger Cat states of dif- 



ferent flow states in an Atom SQUID. These states are 
analogous to the Schrodinger Cat states of different flux 
states demonstrated in SQUIDS [23 EH], and now used 
as qubits in quantum information processors |29| . There 
has been considerable discussion regarding the true num- 
ber of entangled particles in these systems [301 [31] and so 
it would be intriguing to duplicate the experiment with 
a BEC where, in contrast to the SQUID, the microscopic 
physics is almost completely understood. 
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